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We show that filamentous Atmospheric Rivers (ARs) over the Northern Atlantic Ocean are closely linked 
to attracting Lagrangian Coherent Structures (LCSs) in the large scale wind field. LCSs represent lines of 
attraction in the evolving flow with a significant impact on all passive tracers. Using Finite-Time Lyapunov 
Exponents (FTLE), we extract LCSs from a two-dimensional flow derived from water vapor flux of atmospheric 
reanalysis data and compare them to the three-dimensional LCS obtained from the wind flow. We correlate 
the typical filamentous water vapor patterns of ARs with LCSs and find that LCSs bound the filaments on 
the back side. Passive advective transport of water vapor from tropical latitudes is potentially possible. 

PACS numbers: 42.27.De,05.60.-k,92.60.-e,92.60.J- 


Atmospheric Rivers (ARs) transport moisture 
from the Tropics towards northern latitudes and 
they are responsible for extreme precipitation 
and flood events as they approach coastal areas. 
Advection of Lagrangian tracers by geophysical 
flows has shown the presence of persisting trans¬ 
port patterns of Lagrangian Coherent Structures 
(LCS). In this paper, we argue that under some 
circumstances, ARs can be considered attract¬ 
ing LCS in front of which moisture accumulates, 
while the rivers propagate northwards. This be¬ 
havior has been observed both in three dimen¬ 
sional simulations of the wind flow and in terms 
of the two-dimensional vertically averaged water 
vapor field. 


I. INTRODUCTION 


The transport of moisture from the tropics to mid¬ 
latitudes is not continuous and uniform, but rather in¬ 
termittent. More than 90% of poleward water vapor is 
transported by narrow and elongated structures (longer 
than 2000 km and narrower than 1000km), mostly within 
the Warm Conveyor Belt (WCB) ahead of cold fronts 
and within the Low Level Jet (LLJ) of extratropical cy¬ 
clones commonly associated to the polar front (Newell 
et |Zhu and Newell[p"994[ [Ralph and Dettinger[ 


2011| ). These structures, referred to as Tropospheric or 
Atmospheric Rivers (ARs), are defined as elongated re¬ 
gions of Integrated Water Vapor column (IWV) over 2 cm 
and winds stronger than 12 m/s, that transport m oisture 
in the lower troposphere close to the 850 hPa level (Ralph 
[and Dettinger[ |2Q11| [Neiman et al\ |2QQ8[ [Leung and 
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Qian 2009). The advection of moisture by ARs is a key 


process for the Earth’s sensible and latent heat redistri¬ 
bution and has a strong impact on the water cycle of the 
mid-latitudes by generating extreme precipitation events. 
The connection between extreme precipitation and flood 


events has been shown over the Western US Coast (Le- 
ung and Qian[ [2009| [Stohl, Forster, and H.Sodemann[ 


2008) and over Europe (Lavers and Villarini 2013). 


In large scale geophysical flows, advection is often the 
primary process shaping tracer patterns, and effects such 
as turbulent diffusion or other sinks and sources of the 
tracer are secondary. Such tracer patterns have been 
observed in the atmosphere and in the ocean; example s 
include volcanic ash clouds ([Gudmnndsson et al.\ 2012 b 


plankton blooms ( Toner et al ., 2003| Lehahn et al. 


Huhn et al.\ [2012 ), and oil spills (Olascoaga and 


2012). If advection is dominant, Lagrangian Coherent 


2007 


aller 


Structures (LCS) are the relevant finite-time structures 
in the time-dependent flow that determine the deforma¬ 
tion of the fluid and hence the evolution of any advective 
tracer field. The concept of LCSs has its origin in dynam¬ 
ical systems theory. It has been introduced by [Haller and 


Yuan[(2000) and is still being devel oped further ( [Faraz- 


mand, Blazevski, and Hallerj 2014). Attracting hyper¬ 


bolic LCSs are lines evolving with the flow that maxi¬ 
mally attract fluid. In the vicinity of attracting LCSs, the 
fluid is stretched in one direction and compressed in the 
orthogonal direction. Therefore attracting LCSs are the 
cores of filamentous tracer patterns. A standard way of 
estimating the position of attractin g LCSs is the Finite- 
time Lyapunov Exponent (FTLE) ( [Haller 2002 Shad- 
den, Lekien, and Marsden| 2005[ ). The FTLE measures 


the maximum stretching rate among trajectories over a 
fixed time interval derived from neighbor fluid particles. 
Ridges of maximum FTLE val ues are a rob us t method 
to estimate hyperboli c LCS s ( [Haller , 2002 [Shadden/ 


Lekien, and Marsden[ 2005), although some drawbacks 


have been reported (Haller, 2011). 


FTLE and similar versions of the Lyapunov exponent 
have been used to estimate LCSs in different models of at- 





































































































































2 


mospheric flows, such as a zonal stratospheric jet (Heron- 

Vera et al] 2010|), a jet-stream ( 

Tang et al. 2010[), a 

hurricane ([Rutherford et al. 2010 

), transient baroclinic 

eddies (von Hardenberg and Lunkeit 2002) and the polar 

vortex (Koh and Legras 2002 

)■ 


Given that ARs over the Atlantic and Pacific Ocean 
appear as coherent filaments of water vapor with a per¬ 
sistence time of several days up to a week, and given that 
LCSs have turned out to explain the formation of similar 
tracer patterns in geophysical flows, this paper addresses 
the questions: Are ARs associated to attracting LCSs in 
the large scale tropospheric flow? To what extent do ARs 
consist of humid air advected from tropic latitudes? 

Since ARs are defined and studied in terms of two- 
dimensional maps of IWV, we use a two-dimensional 
height-averaged tropospheric flow to detect the main 
LCSs. We focus on selected ARs events that prop¬ 
agate over the Northern Atlantic Ocean and hit the 
Iberian Peninsula. Wind fields and water vapor flux 
fields are retrieved from European Center for Medium- 
Range Weather Forecast reanalysis, ERA-Interim. We 
focus on spatiotemporal patterns in the two-dimensional 
water vapor field formed during an AR event. These pat¬ 
terns are compared with advective structures in the wind 
field, namely LCSs. As mentioned earlier, ARs are th ree- 
dimensional phenomena (Sodemann and Stohl 2013), so 


we have also considered a three-dimensional case study 
to address the limitations of a two-dimensional approach 


(Stohl 1998). 


We find that ARs with a sharp filamentous shape are 
associated to a strong attracting LCS in the wind field. 
This suggests a dominance of horizontal advection in the 
formation process of the filamentous pattern. In these 
cases, the LCS typically marks the northern boundary of 
the AR - a line of high water vapor gradient where humid 
and dry air converge. Filamentous LCSs appear mainly 
in winter, while wide, less coherent ARs structures with¬ 
out a clear attracting LCS develop in summer. 


II. METHODS 

Atmospheric Rivers have been analyzed in terms of a 
two-dimensional vertically integrated flow based on data 
retrieved from the European Center for Me dium-Range 
Weather Forecast reanalysis, ERA-Interim (Dee et al 

2QTT]), 


where Q is the vertically integrated water vapor, $ = 
^ 0 ) are the eastern/northern water vapor flux, g is 
the acceleration of gravity, and 77 is a hybrid vertical coor¬ 
dinate ( [Kasahai^ |1974[ [Simmons and Burridge 1981). 


This coordinate uses the mean sea level as a bottom ref¬ 
erence level and p is the pressure level in the g coordinate. 
u and V are the eastward and northward wind component 
and q is the specific humidity. Data are available with 
a spatial resolution of 0.7° and temporal resolution of 6 
hours. The eastward and northward averaged velocities 
in a longitude-latitude spherical coordinate system ((^, 0) 
are calculated as 


Vf = < X{ip,9,t) >,< > 

.Q’ g 


( 4 ) 


When the velocity fields are vertically averaged in this 
way, the pressure levels are weighted according to their 
water vapor content, such that Vf is representative for 
the dynamics of ARs. 

Along this paper, we will compare the results obtained 
with the previous vector field with those obtained from 
the wind velocity fields at individual pressure levels. 


Vw = [u{ip,9,t),v{tp,9,t),w{ip,9,t)] 


( 5 ) 


Hence, two types of Lagrangian simulations have been 
performed. Either particle trajectories are computed at 
constant pressure levels (2D simulation), or a full 3D sim¬ 
ulation is done. For the last case, the top boundary of the 
domain is located at 300 hPa, while the bottom bound¬ 
ary is at 1000 hPa. Particle trajectories are computed 
by integrating (|^ or ([^ using a 4-th order Runge-Kutta 
scheme with a fixed time step of At = 2 hours, and a 
multilinear interpolation in time and space. A fine grid 
of particles with an initial separation of 1/5° in each pres¬ 
sure level is advected to obtain fields of the Lagrangian 
quantities introduced below. For the two-dimensional 
simulations re = 0 in Eq. <§■ The finite integration time 
is chosen to be r = 120 hours, a typical time scale for 
the formation and propagation of the ARs. 

In order to detect attracting Lagrangian coherent 
structures, we use the Finite-Time Lyapunov Exponent 
(FTLE) field cr(x) which is a measure of stretching about 
fluid trajectories. It is comput ed from the trajectories of 
Lagrangian tracers in the flow (Peacock and Dabiri 2010) 
as 



a{ro,to,T) = — logv^ 

max (C(ro)), (6) 

rl 

where Pmax is the maximum eigenvalue of the right 
Cauchy-Green deformation tensor C = F^F. F is de¬ 
fined by F(ro) = V(r(to + and r(to + r) is the final 
position of the tracers. The time variation of the FTLE 
field has been computed by following the same steps ex¬ 
plained previously, but varying the initial time to in fixed 
steps Ato = 6 hours in order to release a new initial tracer 
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grid. FTLE are computed in forward and backward time 
directions. 

Finite-Time Lyapunov Exponents have been used ex¬ 
tensively to quantify mixing and especially to extract 
persisting transport patterns in the flow, referred to as 
Lagrangian Coherent Structures (LCS). The FTLE at a 
given location measures the maximum stretching rate of 
an infinitesimal fluid parcel over the interval [to, to + r] 
starting at the point tq at time to- Ridges of the FTLE 
field are used to estimate finite time invariant manifolds 
in the flow that separate dynamically different regions. 
Repelling (attracting) LCSs for r > 0 (r < 0) can be 
thought of as finite-time generalizations of the stable (un¬ 
stable) manifolds of the system. These structures govern 
the stretching and folding mechanism that control flow 
mixing. LCSs as FTLE ridges are extracted using the 


Hessian matrix and th e gradient of (j(ro, tp, r) (Shadden, 


Lekien, and Marsden 2005 Sadlo and Peikert 2007). 


To discard weak structures, a threshold criterion associ¬ 
ated to the 90th percentile is used for the FTLE values 
and for the eigenvalues of the Hessian matrix, a fixed 
threshold is used for the minimal connected ridge area. 
For the 3D simulations, the LCSs are estimated using a 
high isosurface of the FTLE field. 


III. RESULTS 
A. 2D results 

Atmospheric Rivers are observed as filamentous struc¬ 
tures in fields of the integrated water vapor column Q, 
Eq. Q. Figure [^a) shows an example of an AR extend¬ 
ing north-westwards from ^ (20°N, 75°W) and trans¬ 
porting water vapor over the Atlantic Ocean towards the 
Iberian Peninsula. The inset shows the daily precipita¬ 
tion rates well above the mean values for this area that 
are caused by this event. The long ridge of large val¬ 
ues of Q connects the tropics with the North Atlantic 
Ocean and represents a continuous AR structure. The 
backward FTLE field is shown in Fig. [^b). The FTLE 
ridge reaching from the Iberian Peninsula to the Gulf 
of Mexico shows a clear relationship to the Atmospheric 
River depicted previously. From a Lagrangian point of 
view, the AR can be considered as an attracting LCS or 
unstable manifold of the flow dynamical system Vf 

Similar results have been obtained for different dates. 
Lagrangian Coherent Structures calculated from the 
backward FTLE fields are compared to the Q field 0 in 
Fig. 0 for four ARs events. Note that regions of max¬ 
imum Q values match LCSs locations. The good agree¬ 
ment of LCSs computed from the flow field Vf with pat¬ 
terns in the tracer field Q indicates that Vf captures the 
main dynamic processes shaping the water vapor field. 

In three cases. Figs. Q(a) to 0 ( c), LCSs and Q 
ridges show a filamentous structure, while for the event 
in Fig. ([^(d), LCSs have a scattered shape. Filamen¬ 
tous ARs have one or two narrow jets growing from low 


latitudes transporting water vapor to medium latitudes. 
However, potential ARs situations with a scattered shape 
do not correspond to well defined water vapor jets. Large 
areas of water vapor surround the jet. The last situa¬ 
tion is typical of summer cases, while the first one occurs 
mainly in winter or beginning of spring. For the winter 
ARs cases, the threshold criteria used to extract LCSs 
from the FTLE fields is insensitive to variations. How¬ 
ever, small variations in the thresholds for the summer 
cases lead to very different LCSs patterns, an indication 
that the LCSs patterns associated to ARs are more co¬ 
herent in winter than in summer. 

In all analyzed winter ARs events, LCSs are located 
behind the river (in the direction of propagation) and do 
not coincide with it. Here the shape of the river is de¬ 
fined as the ridge of the Q field. From a Lagrangian point 
of view, the attracting LCS accumulates water vapor in 
front of the pattern moving towards the east, and a cer¬ 
tain gap in between both maxima, LCSs and Q field, is 
observed. 

These results are obtained using the vertically aver¬ 
aged flow Vf where the averaging is weighted with wa¬ 
ter vapor content. Another natural choice is to define 
a two-dimensional flow on isobaric levels. We compute 
the FTLE fields at different individual isobars using the 
wind Vw 0 and neglecting the vertical wind component, 
w = 0. Figureshows an example of the FTLE fields ob¬ 
tained for the pressure levels 850 hPa and 1000 hPa. At 
850 hPa, ridges of the FTLE field show a close relation¬ 
ship with those shown in Fig. [^b). For other pressure 
levels, we find that the main LCS preserve its continu¬ 
ity, but loses the smoothness of the curvature, showing a 
different shape. 

However, the AR has a vertical extent and it is 
not clear which isobar to choose, such that the two- 
dimensional flow best represents the dynamics of the Q 
tracer, i.e., such that is is most similar to Vf . In order to 
compare trajectories for both two-dimensional flow fields, 
Vf and Vw, we introduce a similarity measure, the mean 
distance error 6. It measures the mean arc length dis¬ 
tance between the final position of the tracers obtained 
for both flows. 


<5 = iV-i^|r, 


- Tv 


( 7 ) 


The subindex of r accounts for the used flow field, and N 
is the number of passive tracers released. The mean dis¬ 
tance error 6 is then calculated for 20 pressure levels from 
300 hPa to 1000 hPa for different ARs events, as shown 
in Fig. A minimum of S is observed for a pressure level 
near 850 hPa for all ARs events analyzed. The minimum 
suggests that the flow field at 850 hPa is closest to the 
water vapor flow Vf. It is consistent with observations 
that the core of ARs is typically located at that height 


( Ralph and Dettinger 2011 Neiman et al. 2008 Leung 
and Qian[ ~ 2009). 
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Figure 1. (a) Atmospheric River event on 15 October 1987 in terms of the integrated water vapor column, Eq. 0 - The ^ flow 
is shown superimposed. The inset shows the accumulated precipitation rates over the Iberian Peninsula [mm], (b) Backward 
FTLE field for the same date and the flow given by Eq. Q 
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Figure 2. Lagrangian Coherent Structures (blue dots) and integrated vapor column Q for four different events. Note the 
coincidence between LCSs and the maximal Q values in the first three cases (winter and spring). 
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(b) March 21,2001 
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(c) April 19,2000 
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B. 3D results 

So far, LCSs have been calculated using the two- 
dimensional flow given by Eq. However, including 
the vertical wind component might significantly alter the 
LCSs obtained from a purely two-dimensional analysis. 
To verify that, we compute the FTLE backward fields 
from a full 3D simulation using the wind field Vw <§■ 
Figure shows the three-dimensional LCSs obtained us¬ 


ing a fixed isosurface from the FTLE field for a winter 
AR. Note the ribbon extending from the Tropics towards 
the Iberian Peninsula. The surface of the LCS is contin¬ 
uous, independent of the pressure level, although slightly 
changes with increasing height. Note that the 3D FTLE 
isosurface follows the Atmospheric River depicted below 
in the figure. This fact suggests that the AR dynamics 
can be considered two-dimensional and the vector field 
defined by Eq. Q accounts for the global dynamics of 
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Figure 3. FTLE backward fields computed using the wind 
field at two different pressure levels, 850 hPa (a) and 1000 
hPa (b), for the ARs event on October 23, 1987. 



Figure 4. Mean distance error defined by Eq. 0 for three 
ARs events at different pressure layers. December 1994 (red 
crosses), August 1988 (blue dots), and November 1995 (green 
triangles). 


the ARs. The vertical FTLE slice shows the presence 
of the atmospheric river at approximately 40°. For this 
latitude, Q diminishes drastically, as the water vapor ac¬ 
cumulates in front of the river as it displaces towards east. 
Note the accumulation of numerous filaments with large 
FTLE values behind the AR due to flow mixing probably 
induced by the interaction between the jet stream and 
the AR (Sodemann and Stohl 2013). 

Even if the horizontal shape of an AR turns out to 
be dominated by the horizontal flow, vertical motions in 
the atmosphere are needed to account for precipitation 
and evaporation processes. These processes may be rep¬ 
resented by the vertical flux of precipitable water on a 
given level. 


(pz = pw {Qv + qi + Qs) ( 8 ) 

where qi and are the amount of water in the 
gaseous, liquid, and solid states, respectively, retrieved 
from the ERA-Interim database, p is the density of the 
air and w is the vertical component of wind. Figure 
shows the vertical flux of precipitable water for two days 
where ARs are present. The vertical flux ^ reaches 
maximum values all along the main attractive LCS bar¬ 
rier, but rather directly in front of the LCS than on top 
of it. This location of maximal vertical moisture flux 
corresponds to ridge of maximal water vapor content, cf. 


Fig.[^b,c. The 3D LCS seems to act as a lateral barrier 
for the vertical transformation processes of water vapor 
that develop in front of it as the AR moves eastward, 
showing that the LCS is a dynamical barrier to precip¬ 
itable water. 

One of the issues that emerges from these findings is 
that high vertical flux also means that Q is not a con¬ 
served tracer, and if we want to characterize vertical mo¬ 
tions of precipitable water active tracers are needed. But 
an implication of this is the possibility that the LCSs 
can identify regions which are not sensitive to the use 
of active tracers for horizontal movements of precipitable 
water. 


C. Origin of air 

Since we have observed stronger attracting LCSs for 
filamentous winter ARs than for scattered summer ARs, 
we also investigate the origin of air parcels in both cases. 
We quantify the displacement of passive particles for two 
ARs, once the AR hits the Iberian Peninsula. To that 
end, an initial grid of particles is advected backward in 
time for r = 5 days. Figure shows the origin of the 
particles in latitude for a filamentous (a) and a scattered 
(b) river. Note in panel (a) that particles coming from 
low latitudes coincide with high values of the integrated 
water vapor column. This potentially gives the possibil¬ 
ity that moisture from the Tropics is passively advected 
northwards, eventually forming the AR. The contribution 
of evaporation and precipitation can, however, not be as¬ 
sessed with this simple simulation of advection. Both 
processes tend to represent significant sinks and sources 
in the water vapor budget. For the unstructured river 
(b), the contours of the original latitude do not enclose 
the ARs structure. Only a very thin zone from low lati¬ 
tudes with nearly zero area is coupled to the AR event. In 
this case, two-dimensional passive advection of moisture 
from low latitudes as a main source can be excluded. 


IV. CONCLUSIONS 

The propagation of Atmospheric Rivers over the North 
Atlantic Ocean, finally hitting the Iberian Peninsula, has 
been studied in terms of Lagrangian tools. Based on 
an integrated water vapor flux obtained from the ERA- 
Interim database, AR events have been identified with 
attracting LCSs, both in 3D and 2D simulations. Two 
different AR events have been characterized. 

On the one hand, narrow filamentous ARs with a fast 
and persistent eastward transport typically develop in 
winter. An attracting LCSs with the same shape as the 
ARs exists in the flow. It acts as a lateral boundary to 
the ARs. This boundary is more clearly visible for the 
3D simulations where the AR is identified as a vertical 
ribbon displacing eastward. As an advection experiment 
shows, for this kind of ARs, the air originates from low 
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Figure 5. Backward three-dimensional FTLE for the Atmospheric River event show in Fig. 0. The isosurface of the FTLE 
field of 0.04 h“^ is shown in grayscale shades while the Q field corresponding to this AR is depicted below. A vertical cut at 304° 
longitude shows the presence of the AR at 40° latitude. Eor comparison, Q values at this longitude are shown superimposed 
over the ETLE field. 

ance of these ARs must be dominated by local sources, 
since the passive transport of moisture from low latitudes 
is practically excluded. 

Our Lagrangian analysis assumes that water vapor 
behaves as a passive tracer. This assumption works 
to assess the deformation in the flow field and obtain 
LCSs that shape tracer patterns in the water vapor field, 
namely ARs. For a quantitative analysis of the water 
vapor balance in ARs, inertial tracers which may ex¬ 
change matter and thermodynamic properties with the 
surrounding flow should be considered. 

Finally, we conclude that the close connection between 
attracting LCSs and ARs should be taken into account 
for future studies and may help to characterize this kind 
of event. Most probably, our results are not limited to 
the Atlantic Ocean, but can be extended to ARs over 
the Pacific Ocean. An extensive analysis with more ARs 
events, including other regions, could help to set up a 
definition of ARs in terms of Lagrangian analysis. 
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(a) March 21,2001 


(b) April 19,2000 



Eigure 6. Vertical flux of precipitable water at 850 hPa for 
two Atmospheric Rivers shown in Eig. §. Eor comparison 
the LCSs detected in both cases are shown as black lines. 



Eigure 7. Contours of latitudes of origin of passive tracers 
advected backward in time for two ARs events. The red, 
green and blue contours enclose particles coming from 25° N, 
21°N, and 17°N. Shaded gray images correspond to water 
vapor concentration Q. 


latitudes and passive transport of water vapor from the 
Tropics northwards is potentially possible. 

On the other hand, for unstructured or scattered rivers, 
mostly occurring during the summer season, a clear at¬ 
tracting LCS seems to be absent. The water vapor bal- 




























7 


with Drs. S. Brands and G. Miguez. 


REFERENCES 

Beron-Vera, F. J., Olascoaga, M. J., Brown, M. G., Kocak, H., and 
Rypina, I. L, “Invariant-tori-like lagrangian coherent structures 
in geophysical flows,” Chaos 20, 017514 (2010) 

Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, 
P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., 
Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bid- 
lot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, 
A. J., Haimberger, L., Healy, S. B., Hersbach, H., Holm, E. V., 
Isaksen, L., Kallberg, P., Koehler, M., Matricardi, M., McNally, 
A. P., Monge-Sanz, B. M., Morcrette, J. J., Park, B. K., Peubey, 
C., de Rosnay, P., Tavolato, C., Thepaut, J. N., and Vitart, F., 
“The era-interim reanalysis: conflguration and performance of 
the data assimilation system,” Quart. J. Roy. Meteor. Soc. 137, 
553-597 (2011) 

Farazmand, M., Blazevski, D., and Haller, G., “Shearless transport 
barriers in unsteady two-dimensional flows and maps,” Physica 
D 278-279, 44-57 (2014). 

Gudmundsson, M. T., Thordarson, T., Hoskuldsson, A., Larsen, 
G., Bjornsson, H., Prata, F. J., Oddsson, B., Magnusson, E., 
Hognadottir, T., Petersen, G. N., Hayward, C. L., Stevenson, 
J. A., and Jonsdottir, L, “Ash generation and distribution from 
the April-May 2010 eruption of Eyjafjallajokull, Iceland,” Sci. 
Rep. 2, 572 (2012) 

Haller,, “A variational theory of hyperbolic Lagrangian Coherent 
Structures,” Physica D 240, 574-598 (2011). 

Haller, G., “Lagrangian coherent structures from approximate ve¬ 
locity data,” Phys. Fluids 14(6), 1851-1861 (2002). 

Haller, G. and Yuan, G., “Lagrangian coherent structures and mix¬ 
ing in two-dimensional turbulence,” Physica D 147, 352-370 
( 2000 ). 

von Hardenberg, J. Fraedrich, K. and Lunkeit, F. Provenzale, A., 
“Transient chaotic mixing during a baroclinic life cycle,” Chaos 
10, 1054-1500 (2002). 

Huhn, F., von Kameke, A., Perez-Muhuzuri, V., Olascoaga, M., 
and Beron-Vera, F., “The impact of advective transport by the 
South Indian Ocean Countercurrent on the Madagascar plank¬ 
ton bloom,” Geophys. Res. Lett. 39, L06602 (2012). 

Kasahara, A., “Various vertical coordinate systems used for nu¬ 
merical wather prediction,” Mon. Weather Rev 102, 509-522 
(1974). 

Koh, T.-Y. and Legras, B., “Hyperbolic lines and the stratospheric 
polar vortex,” Chaos 12, 382-394 (2002). 

Lavers, D. A. and Villarini, G., “The nexus between atmospheric 
rivers and extreme precipitation across Europe,” Geophys. Res. 
Lett. 40, 3259-3264 (2013) 

Lehahn, Y., d’Ovidio, F., Levy, M., and Heifetz, E., “Stirring of the 


northeast Atlantic spring bloom: A Lagrangian analysis based 
on multisatellite data,” J. Geophys. Res. 112, C08005 (2007) 

Leung, L. R. and Qian, Y., “Atmospheric rivers induced heavy 
precipitation and flooding in the western U.S. simulated by the 
WRF regional climate model,” Geophys. Res. Lett. 36 (2009), 
10.1029/2008GL036445 

Neiman, P. J., Ralph, F. M., Wick, G. A., Lundquist, J. D., and 
Dettinger, M. D., “Meteorological Characteristics and Overland 
Precipitation Impacts of Atmospheric Rivers Affecting the West 
Coast of North America Based on Eight Years of SSM/I Satellite 
Observations,” Journal of Hydrometeorology^ J. Hydrometeor 9, 
22-47 (2008) 

Newell, R. E., Newell, N. E., Zhu, Y., and Scott, C., “Tropospheric 
Rivers? - A pilot study,” Geoph. Res. Letters 19, 2401-2404 
(1992). 

Olascoaga, M. J. and Haller, G., “Forecasting sudden changes in en¬ 
vironmental pollution patterns,” PNAS 109, 4738-4743 (2012). 

Peacock, T. and Dabiri, J., “Introduction to focus issue: La¬ 
grangian coherent structures,” Chaos 20, 017501 (2010) 

Ralph, F. M. and Dettinger, M. D., “Storms, ffoods, and the science 
of atmospheric rivers,” Eos, Trans. AGU 92, 265-266 (2011) 

Rutherford, B., Dange, G., Persing, J., Kirby, M., and Mont¬ 
gomery, M., “Lagrangian mixing in an axisymmetric hurricane 
model,” Atmos. Chem. Phys. 10, 6777-6791 (2010). 

Sadlo, F. and Peikert, R., “Efficient Visualization of Lagrangian 
Coherent Structures by Filtered AMR Ridge Extraction,” IEEE 
T. Vis. Comput. Gr. 13(6), 1456-1463 (2007). 

Shadden, S. C., Lekien, F., and Marsden, J. E., “Definition and 
properties of Langrangian Coherent Structures from finite-time 
Lyapunov exponents in two-dimensional aperiodic flows,” Phys¬ 
ica D 212, 271-304 (2005). 

Simmons, A. and Burridge, D., “An Energy and Angular- 
Momentum Conserving Vertical Finite-Difference Scheme and 
Hybrid Vertical Coordinates,” Mon. Wea. Rev 109, 758-766 
(1981). 

Sodemann, H. and Stohl, A., “Moisture origin and meridional 
transport in atmospheric rivers and their association with mul¬ 
tiple cyclones,” Mon. Wea. Rev. 141, 2850-2868 (2013) 

Stohl, A., “Computation, accuracy and applications of trajecto¬ 
ries - a review and bibliography,” Atmos. Environ. 32, 947-966 
(1998) 

Stohl, A., Forster, C., and H.Sodemann,, “Remote Sources of water 
vapor forming precipitation on the Norwegian west coast at 60° 
N- a tale of hurricanes and atmospheric river,” J. Geophys. Res. 
113, D05102 (2008) 

Tang, W., Mathur, M., Haller, G., Hahn, D., and Ruggiero, 
F., “Lagrangian Coherent Structures near a Subtropical Jet 
Stream,” J. Atmos. Sci. 67, 2308-2319 (2010). 

Toner, M., Kirwan Jr., A., Poje, A., Kantha, L., Miiller-Karger, F., 
and Jones, C., “Chlorophyll dispersal by eddy-eddy interactions 
in the Gulf of Mexico,” J. Geophys. Res. 108(C4), 3105 (2003). 

Zhu, Y. and Newell, R. E., “Atmospheric Rivers and Bombs,” Geo¬ 
phys. Res. Lett. 21, 1999-2002 (1994). 


